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ABSTRACT 


Much analysis has been done to di^ on die deformation of helical springs 
due to normal loading. The aim of this study is to design a helical spring diat will 
deform under eccentric loading a desired amount due to a given force. Undo* the 
assun^ons of linear stress strain telationrii4>s, die firing will be designed in 
t«ms of its matoial properties and its geometry. The deformation of the qxing 
wiU be made possible utilizing a Shape Memory Alloy (SMA) active elmient that 
undergoes {diase transformation upon heating above a cotain t e nyerature. Two 
models for spring deformation have been considered. In the first model we study 
the difieiential compression of a spring using SMA wire actuators, and in the 
second model we investigate the bending of an SMA rod placed inside the spring. 


Our efforts were a first step towards the development of a structural skeleton for 


a minimally invasive surgical manipulator. 
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The study of the defornation of helical springs has most comonly been 
limited to those cases due to normal loading IRef. 1, 4, 5, 6, 7, 8] . In 
normal loading, a force is applied through the center of a spring, and the 
displacement of the spring is e3q;>ressed as a function of the load and the 
parameters of the spring involving material properties and geometry. This 
thesis investigates the deformation of helical springs due to eccentric 
loading, with the aim of designing a helical spring that will deform to a 
given shape through the application of known forces. The spring design 
will include the material selection and the selection of helical spring 
geometry in terms of its length, the diameters of the spring and the 
spring coil, emd the number of coils. 

Chapter II provides the backgrovuid needed to study the deformation of 
the two helical spring models, such as the development of Castigliano's 
theorem, and a brief introduction of Shape Memory Alloy's (SMA's). The 
descriptions of the two kinematic models that have been chosen for this 
study are also included in this chapter. Chapter III provides the 
analysis for spring deformation based on the first model, emd Chapter IV 
provides the seune for the second model. Chapter V contains computer 
progreuns that will be used for simulation based design of the proper 
helical spring. 

The motivation of this study is to design a helical spring that will 
be used as the primary structural element of a robotic manipulator for 
minimally invasive surgical applications. In viewing the spring as a 
skeletal element, it is necessary to analyze the deformation of the 
skeleton under the action of external forces, so as to control its 
deformation or motion. The spring deformation will be produced through 
the use of Shape Memory Alloy (SMA) active elements which undergo a phase 





transformation upon haating abova a cartain taiiq;>aratura. Whan thasa SItt*s 
undergo thair phase transformation, they change their shape to a 
predetermined form. During this process they will deform the spring they 
are acting on. By properly selecting the arrangement of the SMA active 
elements about the spring, and controlling their phase transformation, a 
properly designed spring can be made to deform to a desired shape. 
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certain pralininary topics naad to ba raviswsd bafors an analysis of 
tha dafonnation of halical springs dua to accantric loading can ba carriad 
out. This chaptar includas tha daval^pnant of Castigliano's thaoram, an 
introductory axamination of tha propartias of Shapa Manory Alloy's, and 
dascriptions of tha spring nodals v hava chosan for analysis. 

A. CASnOLIAVO'S mORBC 

Castigliano's thaoram is a vary usaful mathamatical tool that can ba 
usad to study tha daformation of alastic bodias undar tha application of 
ganaralizad forcas. Tha dafonnation of tha alastic body is cooputad from 
the strain anargy of tha body. It is an anargy basad approach and tharain 
lies its simplicity. 

Castigliano's thaoram, davalopad by Albarto Castigliano in 1879, is 
a method by vrihich one can determine tha deflection of an clastic body at 
tha point of application of a force. If forcas and F, are axartad on 
an alastic body at two different points, A and B, there are four 
associated deflections [Ref 2] 

( 2 . 1 ) 

4m*<w^s - 2 ) 

(2.3) 

(2.4) 


where the first subscript inplies the point of interest, and the second 





iKqplias tho force of influence. The f'n are constants, and are known as 
influence coefficients. They represent the deflection of one point 
relative to the other. For example, represents the deflection of point 
A relative to point B. These influence coefficients are properties of the 
elastic member. 

Maxwell's law of reciprocity states {Ref. 2, p. 637] 

fai-fsn (2.5) 

which inplies that the deflection at point A due to a unit force applied 
at point B is equal to the deflection at point B due to the application of 
a unit force at A. 

Castigliano developed a method where the deflection due to multiple 
forces acting on an elastic body is obtained as the summation of 
deflections due to forces applied sequentially, one at a time. The final 
result is a set of equations like those found in (2.1) through (2.4). One 
must then superimpose these equations to obtain a series of equations for 
the displacements, or 6's, at the different points where the forces act on 
the elastic body. For two forces F* and F, we obtain 

( 2 . 6 ) 

(2.7) 

Now one needs to determine the work done by each force at each point 
of the elastic body where it acts. The work done by F^^ at point A is 

( 2 - 8 ) 


which, after substituting equation (2.1) becomes 
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(2.9) 


Li]c«wis«, th« work dona by F, at point B is 


( 2 . 10 ) 


If fores Fj^ continues acting on the body while F^ is gradually 
applied, we see that there is additional %«ork done on point A due to F,, 
namely 

( 2 . 11 ) 

Similarly, if force F, continues acting on the body as Fj^ is slowly 
applied, one cem see that there is additional work dons on point B due to 
F^, neunely 

( 2 . 12 ) 

After substituting equations (2.2) and (2.3) into equations (2.11) and 

(2.12), emd by invoking Maxwell's law of reciprocity, equation (2.5), one 
obtains 


(2.13) 

(2.14) 

The total work done on the body, or its total strain energy, if F,^ is 
applied liefore F. is the sunmation of equations (2.9), (2.13), and (2.10): 

lfuFi*2fj^tFg*tggl^) (2.15) 
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If fore* F, is ai^liad first, and than F,^, tha ordar of addition bacaaas 
aquations (2.10), (2.14), and (2.9) 

•NOtai (2.16) 

'nius, tha total work dona on tha body is irraspactive of «Aiathar or F, 
is appliad first. 

Fran Equation (2.16) ona finds that tha displacmiant at point A is 
equal to 

(2.17) 


and the displacement at point B is equal to 


(2.18) 


Castigliano's theorem states that for any force F^ acting on em elastic 
body, the deformation or deflection at the point of application of the 
force Fj is 



(2.19) 


in the direction of F^, where U is the total strain energy of the elastic 
body under the application of forces. 

B. SHAPB imiORT AIiIiOTS 

To cause the deformation of a helical spring one needs to apply eui 
external force. Shape Memory Alloy (SMA) active elements were chosen to 
provide the necessary external forces. An SMA active element has low mass 
and a very high force to mass ratio; this attractive feature allows the 
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miniaturization of tha whoXa structura. A briaf introduction of SMA's is 
iiiq;>arativa to hava a good undarstanding of how tha spring will ba 
daformad. 

A Shape Memory Alloy sMA), is a metallic alloy that is given a 
certain predetermined shape at a high temperature. Once the alloy is 
cooled, it can be deformed, and will remain deformed until heated. Once 
heated above a certain ten^erature, the alloy ‘remembers* its undeformed 
shape 2 md returns to it. 

There are many alloys that exhibit this shape memory effect. Among 
them are Ni-Ti, Ni-Ti-Cu, Cu-Al-Ni-Mn. The alloy is first shaped into its 
desired ‘undeformed* shape at a high temperature, vdien the microstructure 
is in its austenite phase. These ‘undeformed* shapes czm vary greatly, 
but the primary shapes we are considering are those of a thin wire of a 
given length, or a rod with a given circular curve. Once formed, the 
alloy is quenched to allow the microstructure to ccxne into its martensite 
phase. It is now ready to be deformed. 

The alloy can now be deformed by stretching it, bending it, or 
reshaping it by any one of a number of means. It will stay deformed from 
its original shape until it is once again heated up back into the 
austenite region, where it will return to its original shape. It is 
beyond the scope of this thesis to present the microscopic analysis of 
this transformation. It is merely intended to explain what an SMA is and 
how it will be used as an actuator for spring deformation. 

The next two sections provide two spring models that describe the 
positioning of the SMA actuators relative to the spring. The SMA elements 
are in their deformed states initially, emd they revert back to their 
undeformed states once heated above a certain tenperature. During this 
process the spring is deformed. The undeformed shape or the memory shape 
can be given to the alloy by annealing for some time at a fixed 
temperature and then by rapid cooling back to room temperature. In the 
discussion to follow it will be assvuned that the SMA has already been 
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given ics memory shape and accention will be focused on the design of the 
helical spring for achieving the design goal. 


C. MODEX. 1: BELZCAZ. SPRINO tnn>SR tCCamUC COMPIIXSSZOV 

This section considers a spring model consisting of a helical spring 
along with its actuators such that the spring ceui provide two rotational 
degrees of freedom besides a single degree of freed^ for linear 
translation. In this model, shown in Figure 2.1, the helical spring is 
fitted with two end caps. Three SMA wire actuators are attached to the 
end caps just outside of the spring. All of the SMA wires are to be of 
the seune length, spaced 120o apart. During the process in which the SMA 
wires are placed, the spring is given a small initial bias compression. 
This keeps the SMA wires taut and eliminates any slack in the wires. 



.As current is applied to heat one of the SMA's, the active element 
shrinks back to its original ‘undeformed* length. During this process the 
ocher two SMA wires remain in their deformed configuration, and the top 
plane of the spring bends over by virtue of eccentric compression. 












F,, F], emd F, are forces exerted on the spring by the SMA wires, and 
F,, F 4 , and Ft are dumny forces. The deflection of the points of 

application of these forces can be readily obtained using Castigliano's 
theorem. Position vectors, r^'s, from an arbitrary point A to the points 
where the forces are applied are constructed. Angle 6 is a measurement 
taken fr<xn the point where F^ is applied, around in a counter-clockwise 
manner. Using these position vectors, moments and torsions due to forces 
F] through F^ are summed up at A. 

R is the radius of the spring, and L is the length of the spring. E 
in the modulus of elasticity of the spring material, and G is the shear 
modulus of that material. I is the area moment of inertia of the cross 
section of the spring coil, and J is the polar moment of inertia. The 
number of spring coils is n. 

Using these values, the total stain energy of the spring cem be 
calculated. Once this has been done, Castigliano's theorem is invoked, 
and the displacement of the spring at any one of the six points of 
application of the forces can be found. Knowing the relative 

displacements of the different points on the spring coil, the emgle of 
deflection can be confuted. 

For the spring design problem, the angle of deflection is 
predetermined. When a helical spring is chosen, R, I, J, E, G, and n are 
known. From these queuitities the force required to deflect the spring is 
calculated. If this force is one that the SMA wire ceux exert on the 
spring, the spring has been properly designed. If not, some of the 
parameters of the spring, geometry or material, must be chemged and the 
forces recomputed. When the force required to deflect the spring matches 
the force the SMA can exert, the design problem is consisted. 

D. MODKIi 2t HELICAL 8PRIBO UMDER BSNDIMO 

In the second spring model, the assumption is that the spring bends 
under the action of SMA rods. Three SMA rods are attached to the spring 
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incernall/ along ics length and placed 120o apart. Initially the SMA rode 
are in their deformed shape. Vfhen one of the SMA rods is heated it 
regains its undeformed shape and bends the spring in the process. 



Figure 2.2 Helical spring under bending. 


The SMA rod applies a force to each of the spring coils. Each of 
these forces are assumed to have two components, one normal to the coil 
directed towards the center of curvature of the spring denoted F, and the 
other in the tangential direction denoted f. From Figure 2.2, 9 is an 
angular measurement internal to the spring coil measured from the outer 
most point on the spring, moving in a counter-clocicwise direction. R is 
the radius of the spring. The radius of curvature measured to the center 
of the spring is P, and p is the radius of curvature to the inside point 
of each of the coils as they are bent. Angle is the angle of curvature 
of half of the spring. L is the length of the spring, n is the number of 
coils, E is the modulus of elasticity of the spring material, G is the 
shear modulus, I is the area moment of inertia, and J is the polar moment 
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of inertia. For an arbitrary point A on the spring, r denotes the 
position vectors from A to the points of application of the forces, and a 
is the angle subtended by each coil at the center of curvature of the 
spring. 

Consider now only the top half of the spring, since the spring is 
syamtetric emd the bottom half is identical to the top. When an SKA rod is 
bent through the application of current, each coil of the spring is acted 
upon by two forces, F and f. Using these forces, emd the position vectors 
from point A to their point of application, a summation of all of the 
bending moments and torsions at A due to the action of these forces can be 
obtained. The total strain energy due to these moments and torsions is 
calculated, and Castigliano's theorem is applied to compute the 
displacements at each point of application of the forces. 

To reiterate, for the helical spring design problmn, the spring's 
displacement is a given. Knowing the equations for the strain energy, and 
the material and geometry of the spring, one cem work backwards to find 
the force required to bend the spring in this manner. One needs only to 
iterate using the spring material and geometry to find a force 
commensurate with the given displacement. 
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III. DIFORHXTIOH AHALT8IS OW tBK rUtST HILICXL SVRXIIO MODIL 


Racalling th* description of Chs first spring modsi described in the 
preliminaries, one finds that Ft, Pj, and Fj are forces exerted by the SKA 
wires on the spring. F}, F 4 , and F« are duany forces on the spring needed 
to find the displacements at their points of application. Point A is 
arbitrary point on the spring, ^Lnd 9 is the angular measurement from F| 
around in a counter-cloclcwise direction. 



Figure 3.1 Helical spring under eccentric compression. 


There are six position vectors from the six points on the spring where 
the forces are applied to point A. 

/i«Jl(l-co80) i-J?sin05 (3.1) 

fj»-JJ(co30-co86O*)a*Ji(sin6O*-9in0) j (3.2) 


fj » -R (cos6 0* ♦CO80) I*R (sin6 0* -sin0) j 


(3.3) 








(1 l-OlnAj 


(3.4) 


f,*-JI(coa<O*<*«oa0)l-<(aln6O*^slii6)5 ( 3 . 5 ) 

f4«-Jt(cos6-coaC0*)l-A(ainA«alni0*)5 (3.6) 

Sumnation of moeiants about point A yialda 

ir^xPi) -0 (3.7) 


where, 

i-1,2....,« (3.8) 

and where M, and are the x and y coaoponents of the reaction mooients at 
point A. 

After confuting the cross products of (r^ x F^), one cam set the lilce 
vector conponents equal to each other amd find 

(F,♦F 4 ) alnS-i-(F.+F.) (aineo*-alnd) 

(3.9) 

-(Fj-t'Ff) (sin60*'*’SinB)] 


AI^«F[-Fi( 1-C086) ■^F4(l-fCO80) ♦(F,+F,) (CO86-CO860*) 
♦ (Fj+Fj) (cob8-hx>s60*) ] 


(3.10) 


To ccxipute the strain energy, one needs to find the Isending moment amd 
torsion, A and C, on the spring coil. These cam Ise obtained frcxn the 
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reaction noments M, and thrmtgh a coordinate transforaation in the 
following manner: 



(3.11) 


where 


T{9)m 


coflO alnB 

*4lzi0 coa6. 


(3.12) 


Following the coordinate conversions, the moment and torsion equations 
become 


t (-Fi ♦F-4 ♦ (♦Fj-F,) CO860*) slx^ 

♦ (Fa+Fj-F,-Fa)8ln60*cose] 


(3.13) 


t-*[ (Fa+Fa+F,+F4+F,+F,) ♦{-Fa-F,+F,+F*) 8ln60*8inB 
♦ (-Fa+F4+ (-Fa+Fj+Fj-F*) cob 60*) 0080] f^ 


(3.14) 


The shear force at any point of the spring coil can Ise obtained frcai 
the static equilibrium of forces: 

^(Fi+Fa*F,+F4+F,+F,)ic (3.15) 

The expression for strain energy due to shear is 
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(3.16) 



v^«r« Uy is ths snsrgy dus to shsar, f is ths toxm factor (10/9 for a solid 
circla), V is tha shaar forca, L is tha langth of tha spring (2nlt), G is 
tha shaar modulus, A is tha cross-sactional araa of tha spring wira, and 
X is tha intagration variabla danoting tha langth of tha spring coil. 
Aftar making tha siibstitution of J, tha polar mosaant of inart ia, for tha 
circular euraa A fay tha axprassion 



and intagrating with raspact to 0 from 0 to 2nx, tha anargy dua to shaar 
bacomas 




l0tmr*V*R 

1007 


(3.17) 


Ona cem now saa that sinca r is siuch smallar that R, tha r^ tarm will 
dominata tha numarator and maka Uy vary small. Equation (3.17) will be 
CCT^red later with tha strain anargy dua to that of torsion to show this 
difference. 

The strain anargy dua to tha t>anding maamnt is obtained through tha 
equation 



(3.18) 


iidiere U. is the strain anargy dua to tha banding mcxnent, m is the bending 
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BooMnt, L is ths Isngth of ths spring (2nx), «nd BI is ths bonding 
stiffnsss of ths spring. Aftsr intsgration ons obtains 


-2F^r^-FtF^*^*ryF^-P,r,-2FjF,*Fi 


♦#?+F,F4 *Fi) 


(3.19) 


The strain energy due to torsion is obtained through the equation 



(3.20) 


where is the strain energy due to torsion, t is the torsion, L is the 
length of the spring, and JG is the torsional stiffness of the spring. 
Through integration the strain energy due to torsion becomes 

♦5F,F,+3FaF4+2F,F5+3FaF4+3ff ♦5F,F4+3F,F5 (3.21) 
♦2F,F,f3fJ+5F4Fj+3F4F,*3jf+5F,F,+3lf) 


A coii 9 >arison of the strain energy due to shear eund the strain energy 
due to the torsional moment demonstrates how much smaller the strain 
energy due to shear is. Assume that their are eight coils, a spring 
radius of 5 ran, a spring coil radius of 0.5 ran, that forces Fi through F, 
are applied uniformly with a unit magnitude of 1 N, and that stainless 
steel is used as the spring's material. The ratio of strain energy due to 
torsion, U^, to the strain energy due to shear, U,, is 
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(3.22) 


Sine* this ratio is so largs, ons can as# that strain snargy Aia to shaar 
is vary much smallar than tha strain anargy dua to torsion. Whan tha 
spring daforms by banding, tha sasui can ba shoun for tha ratio of tha 
strain anargy dua to banding to tha strain anargy dua to shaar. This 
inqplias that tha strain anargy dua to shaar can ba naglactad. 

Tha total strain anaxgy is simply tha addition of aquations (3.19) euid 
(3.21). 

-F^F^-2F^i-2FiFt*^*ryF^-FyFi-2Firt 
*Fj!*F^F, -F^Ff *Fi*F,Ft ] 

+SFjF, +3 ♦SFjF, ♦3F,F4*2F,F,+3F,F, 
♦3J»J+5F,F4+3F,Fj+2F,F,+3Fi 

♦5F4F5+3F4F4+3^+5F5F4+3Fi1 


By invoicing Castigliemo's thaoram tha six displacements at the points 
of application of the forces are obtained as: 
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(3.24) 

t51^i*6Ji*5F,OF4>2^,+3J>il 

(3.25) 

®>**W I-^i*^a*2F,*F*-/»i-2P-,] 

l31^>5#i*6F,^5/V*3F,+2#il 

(3.26) 

l-2A-i'a*^j*2F4+F,-f’,] 

[2J\*3F,*5F,*6F.*5F,*3F,] 

(3.27) 

•‘"■W f-^a‘2F,-F,*F**2F,*F,] 

*'TW 1^^i*2F,*3Fj*5F^*6F,*5F,] 

(3.28) 
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(3.29) 




As statsd bsfors, F,. F^, and F( ara diaoBy forcas, vihich naans that 
thay do not raally axist, thay ara usad naraly to coa^ta tha 
displacanants at thosa points. Satting Fj - P 4 « F, s 0 ona obtains: 

t6Fj43F,+3F,] (3.30) 


«a- 


4.^ I5Fi45F,^2F,1 


(3.31) 


«3 


onJR* 

2BI 


(-Fi+ 2 F,-F,] 


2JG 


13 ^ 1 + 6 ^,+ 3^,1 


(3.32) 


-2^ 1-2Fi+F 3+F,] *^2^ [2Fi*5F,+5F,l 


(3.33) 


«5 


a “5*1 [-Fj-F3+2F5] [3F3+3F,*6Fj] 


2BI 


(3.34) 


‘•’"W ^*i"2*3**5l ♦■W I5F3+2F3+5F3] 


2 tn; 


(3.35) 


Since all of the SMA actuators will Ije identical, the force they will 
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•x«rt is sxpsctad to b« th« saiM. Consequently, if all of the three 
actuators are activated simultaneously, the deformation at the different 
points on the spring will be the lame. This implies that the spring will 
undergo compression. The interest here lies in the bending of the cop 
plane of the spring. Therefore, all of the three SHh actuators will not 
be activated simultaneously. Now define ^ as the angle by which the top 
plane of the spring bends through when one or two of the SMh wires are 
activated simultaneously. 



The two possible cases to cause the bending of the top pleme of the 
spring are the application of one force, and the application of two 
forces. In the first case one force, F|, is applied. This single force 
will cause displacements to occur at all six points on the spring. Of 
interest are the deflections at points 1 and 4 to determine the bending 
angle of the top plane of the spring. The deflections at points 2, 3, 5, 
and o will be such that the bending will not take place at in any other 
direction. By knowing the displacements at points 1 and 4, one can 
calculate p. With Fi applied one finds 
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(3.36) 




(3.37) 


From Figura 3.3 one csui saa thatsinpla gaonatry ona can saa that 



tanP 



(3.38) 


which means 




(3.39) 


When any cne force, say F^, is applied on the spring, the angle ^ is 


given by 













(3.40) 


The other in^rtemt case is when two forces are applied 
simultaneously, say F, «md F,. In this case the inportwt displacements 
are at points 2 and 5, where 

(3.42) 


This in^lies that for the application of the two forces Fj and Fj 

(3.43) 


Since each of the SMA wire actuators will be identical, the f~rce that 
each will exert on the spring will be equal in magnitude. This means, 
then, if two SMA wires are actuated the resultwt angle of deflection will 
be 




(3.44) 


From these results one finds that if one SMA wire is actuated, or two 
are, the angle of deflection is the same. The difference is the direction 
in which the helical spring bends over. It is now clear that hy 
activating one or two actuators it will be possible to bend the top coil 
of the spring in three separate directions by a positive or a negative 
angle p. The case of actuating all three wires has not been investigated 
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sine* this would producs ths seuna rasult as that of applying a singla 
axial force through the center of the spring, or applying a normal load. 

Now that the effects of helical spring geometry and material 
properties are known in c<Mf>uting the deflection angle p, use of these 
relationships will be instrumental in designing the proper helical spring. 
SMA wires deform by a constant amount of almost 5% of their original 
length in the presence or absence of external forces. However, for this 
property of the SMA to be exhibited repeatedly, the opposing stresses in 
the SMA wires should not exceed a certain value. The design of the 
helical spring will be confuted fr<»> a known value of P that will give the 
force required to produce that deflection. If these forces produce 
stresses that are below the stress limit of the SMA wires, then the design 
is feasible. However, one should try to increase the stresses in the wire 
as high as possible, within limits, such that the spring to be designed is 
not too soft. A spring that is too soft will not be strong enough to act 
as a manipulator for surgical applications. 
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IV. OnrORlATIOH JOaLTSIS or TBM SMCOKD HELICAt SrUMO MOOIL 

The second helical spring model considered is one ^lrt^ere the spring is 
bent over directly by an SMA rod. Three rods will actually be placed 
inside the spring 120o apart to give the spring a £ull range o£ motion. 
Here, consideration is given to only one SMA rod for the development of 
the model. 



Figure 4.1 Helical spring under bending. 


The only major assumption made is that the length of the spring is the 
s£une in the bent over configuration as in the undeformed configuration. 
With this in mind, certain quantities are defined. Let R be the radius of 
the spring and L be the length of the undeformed spring. In the deformed 
configuration, it is desired that che spring bends into the shape of a 
circular arc. Let P be the radius of curvature of the central eucis of the 





spring, and p tha radius of eurvatura of tha innar part of tha spring. 
Lat ^ ba tha angla batwaan tha top coil of tha spring and tha cantral 
coil, ^ ba tha angular changa of spring coils, 0 ba tha angular 
maasuranant along tha langth of tha coils of tha spring, a ba tha angla 
subtandad fay two adjacant spring coils, and n ba tha total nusibar of 
coils. Fran thasa dafinitions. ona saas frcn Pigura 4.1 that tha 
following idantitias hold: 



(4.1) 


(4.2) 

n 

p-P>Jf (4-4) 

In this modal thara ara thraa diffarant coordinata framas that ara 
used. Tha traditional t, and Jf frame is fixed on tha top coil of tha 
spring. Tha f', , and Jt' frame is relative to each succeeding coil. 

The e„ and freuna is a normal and tangential tremsfomation of tha f, 
and It frame. 

Suppose that as tha SMA bands over, it exerts two forces on each 
spring coil it touches, F in the normal direction, and f in tha tangential 
direction. The line of action of tha normal forces passes through tha 
center of curvature. For tha ease of coii 9 >utation, it is helpful to 
translate the point of application of the normal force along its line of 
action to the center of curvature of the spring. Point A, at the top 
center of the spring, is the point where position vectors from all other 
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points on ths spring will bs constructsd to. In considsring only ths 
first coil, thsre are two position vectors from each point on that coil, 
one corresponding to the normal force F, and one to the tangential force 
f. These position vectors have been defined this viay so that they will 
apply to each coil for these forces. They are 

(4.5) 

'■ (P-J{ooa6) I '^ItairfbS '-pi (4.6) 

The primes on the i euid j vectors indicate that they are in the spring 
coil's frame of reference, where the unprimed vectors are in the fixed 
frame of the top coil. These vectors will give the position fron any 
point on the first coil relative to point A. After recognizing that 

Fjl ■ F’iCOs4l'’-F,8in^it' (4.7) 

fife » fi8in#l'*f,co84ic' (4.8) 

the total moment about point A caused by these forces cem be calculated )3y 

♦ (li 'xli) (4 • 9) 

For the second coil the position vectors remain the seune, but the 
forces acting on this coil are F^, fj, Fj, euid fj where 

^,*f\l*FiC08 (a+4) 1'-Fi8in (a+4) (4.10) 


£,*fil»fi8in(a+#) 1'♦f,co8 (a-^) £ 


( 4 . 11 ) 








(4.12) 


(4.13) 


Ndw thm total monant about point A is 

J4,- (♦ (fi 'xf,) ♦ (/iX#,) ♦ (fj 'xf,) (4.14) 


By procading in a li)ca mamnar, ona arrivas at tha )c-th coil, vihara 
this coil is aetad upon by tha forcas Fj, Fj, —, F^, and fi, f,, ..., fk. 


The e forcas tales tha fora 

^j-F^cos [ (A-1) f'-F,sin t (*-l) S”^] if' (4.15) 

/,«FaC081 (*-2) «-»41 f'-FaSlnC (it-2) if ' (4.16) 

I 

P,i^Ffpoa^I'~F^ia^' (4.17) 

fi«fi8in ( (*-l) a*4] I'*f^coB ( (A-l) if' (4.18) 

/a»/a8in [ (A-2) a*#] f'•••f,co8 C (*-2) 844 ] if ' (4.19) 

{ 

/fc»faSin#l'+f|^08#ic' (4.20) 


Now the total moment about point A becranes 
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(4.21) 




or, aftor caaputing tha cross products 


iteloA (r^sln I (Jt-i) *44] -^icos [ (Jt-i) ) i' 

*■ [ (P-*JBeoa6) (#iSinI (Jt-i) c-^] 

-f^s I (Jt-i)«441 ) ♦/iP]i' 

♦JlsliiA (F^cos [ (Jt-i) «<»4] *fimin [ (Jt-i) ««4]) J^') 


(4.22) 


Ths sasisst way to worlc with this mcoant is to braa)c it into parts as 
follows 


(-Jlsinisin I (Jr-i) «*41) 
(JZslz^Bcos [ (Jt-i) «44])) 


(4.23) 


"i5 ^ (jf-i) •♦♦1; 

*fi ((P-JIsini) cos [ (Jr-i) «-»4] ~P)) 


(4.24) 


(f^ (-JtslnOcos ((Jr-i) ■♦♦I) 
(-JtsinSsin[ (Jt-i)■'»4])) 


(4.25) 
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A coordinat* transformation is now in ordsr to transform tha meawnts 
from ths coil fixsd frams to a nozmal~tan 0 mtial frams. Ihs 
transformation matrix T(9). from aquation (3.16), will ba usad such that 



coa0 sinftW'V 
-slnA cosaJIk/ 


(4.26) 


Tha raooksnt aquations (4.23) through (4.25) baccem 

m-g {fi (-Phlnasln I (Jr-i) ) 
*fj (PsinAcos ((Ir-i) ««4] -psinA)) 


< (Jl-FcosA) sin I (Jt-i) s-*4]) 
*f^ (- (Jt*^osA) cos [ (Ar-'i) s-»4] *psinA)) 


(Fj (-JhiinAcos I (Jc-i) s*#]) 
*fi (-JZsinAsin [ (Ir-i) s-^])) 


With each of these moments in the proper frame of reference, the 
strain energy due to each of them must be calculated. With the subscript 
Ic on each of the energies to remind one that this is the strain energy of 
the ]c-th coil, the strain energy due to each of these moments is given by 
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(4.30) 





-fil 


■( 


2je 


(4.31) 



(4.32) 


Sine* Mch of thm nooMiit t«nna ar« •xpraasad as a sun of Ic quantitiss, 
ths intagrals involving tha squars of tha nooiants bacooias quita 
cnjnbarsaiaa. To simplify tha caaq>utation, ona naads to usa matrix algabra. 
First dafina tha incramantal forca ^wetor as follows 
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'i;;' 


AFM 


AF^ 

aa 

A/, 


\,At^ 


€■“ 


(4.33) 


Than define 


( >PblnO[(lc>-i)«441 ^ 
>BilaBsl&((Jt»2)c44] 


-Fiixidaln# 


cl* 


(4.34) 



(4.35) 


It can now Ise shown that 
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(4.36) 



Likawis* for th« torsional 


»nt, dsfina 


:«1 


AH 


(Jt-Fcbs6} sin [ (Jr-2) i 


(H^PcosS) sln^ 


(4.37) 


- (Jl-Pcoa6) cos [ (Jc-l) s-»4] -pcos6^ 

- (JI-Pcos6) cos [ lk~2) s«#] -pcos6 


(4.38) 


(JI-Pcos6) cos^'PcosO ) 


Then, 



(4.39) 


For the moment in the z' direction, define 
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Than, 




(4.40) 


-iaiiiKln [ (Jr-2) «■»#] 


-JteidBsln^ 


(4.41) 




(4.42) 


With the moments defined lay Equations (4.36), (4.39), and (4.42), one 
cam obtain their squares as follows: 


HaiiAj) -af 


(4.43) 


- *(B,iE^) -AF 

W 


(4.44) 
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(4.45) 






■A^ 


Now d«fin« 


/jl^\ A]A^ > A}A, A^ I Au 

A^- -(AiiA,)- - - - --€ !**«»* 

W [aA J A^ |a,j I A^ 


one arrives at 


I^-AJ^'-A-A^ 


where the elements of are of the form 

A„ (l,i) -F 8in*« siftt (*-i)«•»♦) sin [ (A-J) 


the elements of Aj 2 are of the form 

Ai,(i,i)—F sin^ 8in[(A-i)s-»A]coB[(A-i)c441 
+Pp8in*e 8ln C (A-i) 844] 


Aai{i,i)-Aja(i,J)>- 


and the elements of A 22 are of the form 


(4.46) 


(4.47) 


(4.48) 


(4.49) 


(4.50) 
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(i.i) ■f^ala^eofl [ GKM ( 

-l^sln*6(coa ((Jr*i)«<»4] *com [ (Jr>i)c«4]) 
♦p*aln*S 


(4.51) 


Similarly, by dafining 



*1^1 fBii « 


BA Ai < B„ 


(4.52) 


cna can aaa that 

I^.AF^’^AF (4.53) 

whara alamants of Bu taka tha form 

a,, (i, • (a-PcoaS) *ain [ (*-i) a*4] ain [ (Jt-J) a-^] (4.54) 

the elements of B,, taka the form 

au(i,^)--(a-Pco80)*eint(Jt-i)a+^]coa[(it-J)a*^] 

(4.55) 

- (a-Fcoa6) pcoaSaln [ (k-i) a^^] 

(4.56) 


and the elements of B^j take the form 




Bja ii.j) ■ (J2-Fco86) *co8 I (Jc-1) a-»41 co8 [ ik-j) 

+ (JZ-Pco 86) pco86 (co8 [ (Jc-i)«♦♦] •k:o 8 t Ik-J) ) (4.57) 

♦P*co8*6 


By defining 



iqq, 

ic^c^ 


Pii ‘ Pu 
^2 * ^2. 


(4.58) 


one can show that 

(4.59) 

where elements of Cji take the form 

C,i (1, J) «ie*8in»eco8 ((Jc-i) cos t (Jr-i) k-^] (4.60) 

the elements of take the form 

C„ U.j) -J?*ain*68in [ (Jt-1) cos I {k~j) a^] (4.61) 

C„(i.J)-C„(i.J)»- (4.62) 

and the elements of C 22 take the form 

Cia (1, J) -J?*8in*e8in ( (k-i)a*^] sin [ (Jc-J) a 44 ] (4.63) 

Now that the squared moments are in contact form, the integrals can 
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b« avaliMtcd for th« caaaputation of tha strain saargias. Tha strain 
anargy in tha k-th coil dua to tha banding aiaaiant 1^ is 



(4.64) 


whara ona must avaluata tha intagral of A tarn by tarn. Ibis maans that 
if ona dafinas 


2 « 

* 2§r/^ (ic-i) Bin ( (k-J) «+♦] dB 

0 

'• 1 •I'O* t* >1 


. sinaain [a Uk-i-J*!) ] 
1 2SI/ a («*-4«*) 

--^sinCaCi-^)] 










and 

■(-=—=){-•ln»cos [« (2k-i~J*l) 1 

\ 2Bll a («*-4ii*) 

♦■2^coa [c U-J) ] ♦«p* 

♦_i452fe_8in|i\co* [i (2Jt-i-J+l) 1 COB [-1 (i-i) 11 

«(«*-16ll*) \2/ 2 2 


then, the energy in the k-th coil due to the nx»nent Mn is 

~ -AF (4.69) 

kr * -*»»: 

The energy in the k-th coil due to the torsional moment takes the 

form 


(4.67) 


(4.68) 
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(4.70) 


wh«r« a tern by tarm intagration, as bafora, yialds Equation (4.75), whara 
tha alamants of Bu'> Bji'. and B 33 ' ara dafinad as follows: 


42St- [•(2Jc-i-^*l) 1) 
^ a «a-«a «(aa-4s*) } 


(4.71) 


(4.72) 

V « «*-«* «(«*-4s») ) 

J 4 n (li»>e«»)^ . j j, 

\ a(«*-l6is») «»-4«*J ^2/ * 


(4-73) 

and 
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Thus the strain energy o£ the k-th coil due to the torsional manent is 

Ai' ‘ Bu" 

Ce ~ ~ ^AF (4.75) 


The strain energy o£ the k-th coil due to the moment M,. is 



(4.76) 


where the elements of the submatricies Cjj', Cj,', Cj,', emd Cj,' are defined 
in the following nuuiner: 
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4-Hcos [«(i~i)l) 


q.-U.J) 1 


and 

♦|cOBl«(i-i)]l 


Thus the energy of the k-th coil due to the moment M,. is 







(4.81) 




Ax' ‘ Ax 
Ax' * Ax 



Th« total anargy of tha k-th coll of tha halical aprlng is tha 
svmnation of tha anargias dua to tha thraa noaiants, or 

A'A k*A. k*A\ s-AF'*U -A^ (4.82) 


Zf ona dafinas 


DAA'*B'*€r 


(4.83) 


tha incramantal anargy of tha halical spring up through tha k-th coil now 
t>ecanes 



(AF***X)*AF)j 


(4.84) 


Castigliano's thaorem is now invokad to detamine the incronental 
displacesiants of aach of tha coils in tha normal auid tangantial 
directions. 
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.. *2Djkr^ 


(4.tS) 



<2DuAFx*I^rt*2Dgjtry*... *2DgJkF^ 


(4.86) 




(4.87) 


which simplifies to 


N 




1 


W 



fl„ 2D„ 2D„ • • • 2D^ 
2^1 OgM 2^ * • 2lJ 

i i ! i 

2D^2D^2D^ • • Dtt, 


AF, 

i 

AF^ 


(4.88) 


Now define 


A4^ 


Ou 

2Du 

2i)„ • 

• • 2D,» 

2D„ 

D„ 

2D„ • 

♦ • 2D^ 

i 

i 

i 

i 


2Am 

21V, • 

• • ^ 


(4.89) 


then the final important relationship IsecoiBes 
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(4.90) 


With this rslationship, knowing ths physical gsookstry and aatarial 
propartias of tha halical spring, and spacifing a dasirad displacasiant, 
ona can datarmina tha forca raquirad by tha SHA activa alMsant to band tha 
spring. Tha matrix A contains a congloataration of inforaiation basad on 
all of thasa known quantitias. Onca thasa givan quantitias hava baan 
spacifiad, it is a sinpla mattar to datarmina tha raquirad forca matrix F. 
If tha forca raquirad cannot ba achiavad by tha SHA, itaration is raquirad 
for tha spring's gaomatry or matarial propartias. A program is anclosad 
that will aid in tha dasign of a halical spring undar thasa conditions. 
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This tlMsis investigated the defontation of helical springs under 
eccentric loading using two different nodels. The goal %«as to develop an 
algoritln for coaputing the deflection of a helical spring under a known 
force. The desired final shape of the spring tots that of a spring bent 
over in a circular fashion. The t%io sodels used achieved this by two 
different methods, both utilizing SKh active elsnents in two different 
configurations. 

The development of this study t#as to obtain relationships betvreen 
force and displacement so that one could predict or control the 
deformation of the spring. A spring under normal loading is governed by 
an equation such as 

(5.1) 

where force F is proportional to displacement S. This proportionality 
consteuit, k, is based on the spring geometry and material properties. In 
the development of this study, similar relationships were obtained, but 
with a more conplex form of proportionality factor. This factor is 
dependent on the spring configuration, and is not a constant. The 
proportionality factor is also a function of the geometry and material 
properties. 

This study has been a first step in the design of a robotic 
manipulator to be conposed of a series of helical springs that utilize SMA 
actuators to control the manipulator. In the course of future research, 
SMA actuator design is needed, emd in assmnbling the helical spring with 
the actuators. Once the actuators are designed, the results acquire by 
this study can be employed to match a spring's geometry and material to 
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th« fore* that can ba ganaratad fay tha SMX to band tha spring in a dasirad 


mannar. 

Whan tha dasign procass is ccnplatad, a halical spring and actuator 
systam will ba raalizad, and can ba inplamantad in a mininally invasiva 
surgical nuuiipulator. Tha conclusions of this study can also ba usad in 
any othar situation whara an accantric load is placad on a halical spring 
causing it to bend over. 
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nooBin A» noM 

Th* first algorithai dslsvppsd %iu for tho first nodol, whsrs ths 
spring was dsformad fay applying an accantrie load axially on tha spring. 
It is a FORTRAN program vAiara matarial pr^artias, spring gaoawtry, and 
dasirad daflaction angla ara raquirad as inputs, and tha forca raquirad to 
dafoxm tha spring is tha output. 

PROGRAM THESIS 

* 

* This progreun is for tha spring dasign usad in tha manipulator skalton 

* for minimally invasiva surgical applications. Givan tha physical 

* dimensions and matarial properties of tha spring, and tha dasirad 

* deflection where tha forca will be applied, it will solve for 

* the deflection of the spring on the opposite side of where tha force 

* is applied, the force required, and for tha angla of daflaction. 

* 

* Inputs: Bending stiffness of coil El, torsional stiffness of the coil 

* JG, number of coils n, radius of the spring R, radius of the spring 

* coil rc, length of spring, emd deflection d where the force will be 

* applied. 

* Outputs: deflections d, deflection angle B, and force F. 

* 

REAL F, E, I, J, G, R, rc, L, pi, c, El, JG, dF, dOPPF, B, Z 
INTEGER n, m, k 
CHARACTER*20 MATL 

PRINT* 

PRINT*,'THIS PROGRAM WILL AID IN THE DESIGN OF A SPRING FOR USE' 
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PRINT*, * IN THE SKBLTOH FOR MINIMALLY INVASIVE SURGICAL 


PRINT*, ' APPLICATIONS.' 

* 

pi s 3.1415926 

* 

5 PRINT* 

PRINT*, 'SELECT AND ENTER THE TYPE OF MATERIAL YOU WISH TO USE' 
PRINT*,' (i.e. STAINLESS STEEL): ' 

READ*,MATL 

PRINT* 

PRINT*,'ENTER MODULUS OF ELASTICITY (E). IN ENGLISH UNITS' 
PRINT*,' (i.«. 28.0«6): ' 

READ*, E 

PRINT*,'ENTER SHEAR MODULUS (G), IN ENGLISH UNITS' 

PRINT*,' (i.e. 10.6e6): ' 

READ*, G 
PRINT* 

PRINT*,'ENTER 'I" IF YOU HAVE DIMENSIONS IN ENGLISH UNITS (in), 
PRINT*, ' OR ■2" IF YOU HAVE DIMENSIONS IN SI UNITS (mn) : ' 
READ*, Ic 

PRINT*,'ENTER RADIUS OF SPRING; ' 

READ*, R 

PRINT*,'ENTER RADIUS OF SPRING COIL; ' 

READ*, rc 

PRINT*,'ENTER LENGTH OF SPRING: ' 

READ*, L 

PRINT*,'ENTER DESIRED DEFLECTION: ' 

READ*, dF 

10 IF (k.EQ.2) THEN 
R = R*(0.001)*39.37 
rc = rc*(0.001)*39.37 
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L . L*(0.001)*39.37 
dP « dP*(0.001) *39.3-’ 

EMDIF 

PRINT* 

PRINT*,'ENTER NOMBER OF SPRING COILS: ' 

READ*, n 
PRINT* 

PRINT*. 'ENTER ‘1* IF (»<E FORCE IS APPLIED,' 
PRINT*,' OR •2‘ IF TWO FORCES ARE APPLIED: ' 
READ*, m 


I = 0.25*pi*(rc**4) 

J = O.S*pi*(rc**4) 
c s n*pi*(R**3) 

El - E*I 
JG = J*G 
IF (m.EQ.l) THEN 
F = dF/(c*(l/(EI) + 3/(JG))) 
dOPPF = c*F*(-l/(EI) + 1/(JG)) 

B = (ATAN((dF - dOPPF)/(2*R)))*(180/pi) 

ELSE 

F = dF/{2*c*(l/(2*EI) + 5/{2*JG))) 
dOPPF = 2*c*F*(-l/(2*EI) + 3/(2*JG)) 

B = (ATAN((dF - dOPPF)/(2*R)))*(180/pi) 

ENDIF 

PRINT* 

PRINT*,'FOR MATERIAL '.MAIL 
PRINT* 

PRINT*,'DEFLECTION AT POINT OF APPLICATION (in); ',dF 
dF = dF/(0.001*39,37) 
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PRINT*, 'DEFLECTION AT POINT OF APPLICATIW (BBi) : '.dF 
PRINT* 

PRINT*,'DEFLECTION AT POINT OPPOSITE APPLICATION (in): ',dOPPF 
dOPPF = dOPPF/(0.001*39.37) 

PRINT*,'DEFLECTION AT POINT OPPOSITE APPLICATION (mn) : ',dOPFF 
PRINT* 

PRINT*,'DEFLECTION ANGLE (degrees): ',B 
PRINT*,'FORCE NECESSARY (lbs): ',F 
PRINT* 

PRINT*,'IF YOU WISH TO BEGIN J^SAIN, TYPE ANY NUMBER AND HIT' 
PRINT*,' RETURN. ELSE ENTER *99* TO QUIT. ' 

READ*, Z 

IF (Z.NE.99) GO TO 5 

END 


The second algorithm developed was for the second model, that of a 
spring bent over directly by an SMA active element. It is a MATLAB 
progreun where one must input the geometry and material properties of the 
helical spring, as well as the final desired deflection angle of the 
spring. It generates a proportionality matrix whose size is based on the 
nvunber of coils of the spring, and computes the force required to bend the 
spring. It plots force versus euigle deflection to obtain a relationship 
between the magnitude of the force required and the amount of deflection. 

% thesis.m 
% 

% This progreun is for the spring design used in the manipulator slceleton 
% for minimally invasive surgical applications. Given the physical 
% dimensions and material properties of the spring, and the desired 
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% deflection of the spring, it will solve for the force required to bend 
% the spring using an SMA rod inside the spring. 

% 

% Inputs: ■ Spring length (L), radius (R), coil radius (r). Modulus of 
% Elasticity (E), Shear Modulus (G), number of coils (n), 

% desired angle of deflection (phi). 

% Outputs: Force required to bend spring (F). 

% 

clear 

% 

L = input('Length of spring in millimeters '); 

R = input('Radius of spring in millimeters '); 
r = input('Radius of spring coil in millimeters '); 

E = input('Modulus of elasticity in GPa '); 

G = input('Shear modulus in GPa '); 

n = input('Number of spring coils (must be even number) '); 
phifinal s input('Final desired deflection angle in degrees '); 

% 

L = L/1000; 

R = R/1000; 
r = r/1000; 

E = E*le9; 

G = G*le9; 

phifinal = phifinal*pi/180; 

% 

I = 0.25*pi*r^4; 

J = 0.5*pi*r''4; 
mom = R/{2*E*I); 
tor = R/(2*J*G); 

% 

for i = l:n 
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initfd) m 0; 
end 

f = initf'; 

% 

count = 0; 
dradius = 0; 
radius(1) = 50; 
for z s 1:200 

anglephi = L/(2*radiua(z)); 

if emglephi <= phifinal, count = count +1; end 
dradius = -0.05*radius(z); 
radius (z-fl) = radius (z) + dradius; 
end 

% 

dP = 0; 

P(l) * 50; 
for q * l:coimt 
% 

for i = l;n 
for j = l;n 
D(i,j) = 0; 
end 
end 

% 

ro s P{q)-R; 

phi(q) = L/(2*P(q)); 

a = (2*phi(q))/n; 

cl = (2*pi-3)/{a*(a^2-4*pi-2)); 

c2 = pi/2; 

C3 = P(q)''2; 

c4 = (8*pi^2)/(a''2-4*pi^2) ; 
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c5 . (64*pi*3)/{a*(«^2-16*pi*2)); 
c6 » pi*(R''2+{P(q)''2/2)); 

c7 = ({R''2*pi)/a - (2*pi*a*R*P(q))/(a^2-pi*2) ... 

+ pi* (a''2-2*pi'‘2) *P(q) *‘21 (a* (a'‘2-4*pi*2))) ; 
c8 » {(4*pi*(a^2-8*pi*2)*P(q)*ro)/(a*(a^2-16*pi*2)) 
- (4*pi*a*R*ro)/(a''2-4*pi''2)); 
c9 = {32*pi''3)/(a*(a'‘2-16*pi'‘2)); 
clO s sin(a); 
cll s sin(a/2); 
cl2 s inc»i*R''2; 
cl3 * cl*c3*cl0; 
cl4 = c2*c3; 
cl5 = P(q)*ro*c9*cll; 
cl6 * c7*cl0; 
cl7 s cl*clO; 
cl8 = P(q)*ro*c5; 

for k = l:n/2 
for j * l:k 
for i = l:k 

cl9 = cos(a*(2*k-i-j+l)); 
c20 = sin(a*(2*k-i-j+l)); 
c21 = cos(a*(i-j)); 
c22 = sin(a*(i-j)); 
c23 = sin(a*(k-i+0.5)) ; 
c24 = sin(a*(k-j+0.5)); 

All(i,j) = inoin*(cl3*cl9 + cl4*c21); 

A12(i,j) = inoin*(cl3*c20 - cl4*c22 - cl5*c23); 
A21(i,j) = mom*(cl3*c20 + cl4*c22 - cl5*c24); 


53 










% 


A22(i,j) 


* mom*(-cl3*cl9 + cl4*c21 ♦ pi*ro^2... 

■*. cl8*cll*co»{(*/2)*(2*k-i-j+l) )*co«( (a/2)*(i-j))) ; 


% 


% 


Bll(i.j) 
B12(i,j) 
B21(i,j) 
B22(i.j) 


s tor*{c6*c21 - el6*cl9); 
a tor*(c6*c22 - cl6*c20 + c8*cll*c23); 
a tor*(-c6*c22 - cl6*c20 + c8*cll*c24); 
a tor*(c6*c21 ♦ cl6*cl9 - c8*cll*(cos(a*(k-i+0.5))... 
+ cos(a*(k-j+0.5))) + pi*ro‘*‘2); 


ClKi, j) 
C12(i.j) 
C21(i.j) 
C22(i,j) 


cl2*(c2*c21-cl7*cl9); 
cl2*(c2*c22-cl7*c20); 
cl2*(-c2*c22-cl7*c20); 
cl2*(c2*c21+cl7*cl9); 


and 

end 


% 


dll a All+Bll+Cll; 
dl2 a A12+B12+C12; 
d21 a A21+B21+C21; 
d22 a A22-«-B22-»-C22; 
for i = l:k 
for j a l;k 

D(i,j) = D(i,j) + dlKi.j); 

D(i+k,j) a D(i+k,j) + d21(i,j); 

D(i,j+k) a D(i,j+k) + dl2(i,j); 

D(i+k,j+k) a D{i+k,j+k) + d22(i,j); 
end 
end 
end 

D a 2.*D; 
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for i « l:n 

D{i.i) * 0.5*D(i.i); 

•nd 

% 

for i s l:n/2 

th«ta(i) » {n/2 - (i-l))*a; 
daltan(i) s (1 - cos(thata(i)))*dP; 
doltat(i) s (sin(thata(i)) - thata(i))*dP; 
end 

% 

delta s [deltan';deltat']; 

F * inv(D)*delta; 
for i = l:n/2 

fx(i) = F(i)*coa(theta(i)) > F(n/2-t-i) *ein(theta(i)); 
fy(i) s F(i) *sin{theta{i)) + F(n/2-i-i} *co8(theta(i)); 
end 

fxy = Ifx';fy'l; 
f * f + fxy; 
force(:,q} = f; 

PHI(q) = phi(q)*(180/pi); 
dP = -0.05*P(q); 

P(q+1) * P(q) + dP; 
end 
% 

clg 

for j = l;n/2 

plot(PHI,force(j,:}) 

title('X-Force vs Radius Angle Phi, for Coil #') 
xledsel (' Phi (degrees) ') 
ylabeK'Force (N) ') 
meta thesplot 
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•nd 


for j « n/2-fl:n 
plot(PHZ,forc«(j,:)) 

titl«('Y-Fore« vs Radius Angls Phi, for Coil #') 
xlabsK'Phi (dsgrsss)') 
ylabsl('Fores (N)') 
meta thssplot 
and 


Ths final sl«Bsnt is a sat of plots using tha pravious program. Thay 
ara forca varsus daflaction angla graphs ganaratad by tha program tising 
tha inputs of a 20 an long spring, with a radius of 5 bb, coil radius of 
0.5 mm, and sight coils. Tha matarial chosan ms stainlass staal with a 
modulus of alasticity of 190 GPa, and shaar modulus of 73 GPa. Tha final 
dasirad daflaction angla ms 30a. Tha sight plots ara for tha t^ half of 
the spring. Tha forces in tha bottom half ara going to ba identical as 
tha spring is synmatrical. Of tha sight forca/daformation plots, four ara 
for normal forces, and four ara for tangential forces. 
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